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Abstract 

A study of Z boson pair production in e + e~ annihilation at center-of-mass energies near 183 GeV and 
189 GeV is reported. Final states containing only leptons, (l + £~£ + £~~ and t r t~vv), quark and lepton 
pairs, (qq£ + £~ , qqvV) and the all-hadronic final state (qqqq) are considered. In all states with at least 
one Z boson decaying hadronically, qq and bb final states are considered separately using lifetime 
and event-shape tags, thereby improving the cross-section measurement. At y/s = 189 GeV the Z-pair 
cross section was measured to be 0.80 ^^(stat.) -ao5( s Y s t-) pb> consistent with the Standard Model 
prediction. At = 183 GeV the 95% C.L. upper limit is 0.55 pb. Limits on anomalous ZZ7 and ZZZ 
couplings are derived. 
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1 Introduction 



The study of the process e + e~ — > ZZ has recently become possible since LEP now operates at center- 
of-mass energies above the threshold for on-shell Z boson pair production. In the Standard Model, 
the process e + e~ — > ZZ occurs via the NC2 diagrams Q shown in Figure |]. The Z-pair cross section 
depends on properties of the Z boson (m^, Tz and the vector and axial vector coupling of the Z 
to electrons, gy and g\) that have been measured with great precision at the Z resonance 0. The 
expected Z-pair cross section increases from about 0.25 pb at \fs = 183 GeV to about 1.0 pb at 
yfs = 200 GeV, but remains more than an order of magnitude smaller than W-pair production. In 
contrast to W-pair production, where tree level WW7 and WWZ couplings are important, no ZZZ 
and ZZ7 couplings are expected in the Standard Model. However, physics beyond the Standard Model 
could lead to effective couplings || which could then be observed as deviations in the measured Z- 
pair cross section from the Standard Model prediction. Such deviations have been proposed in the 
context of Higgs doublet models || and in low scale gravity theories ||. In this paper we report on 
measurements of the NC2 Z-pair cross section, including the extrapolation to final states with one 
or both Z bosons off-shell. These measurements, along with the angular distribution of the observed 
events, are then used to extract limits on possible ZZZ and ZZ7 couplings. 

In Section [2| we describe the data sets used and the Monte Carlo simulation of signal and back- 
ground. In Section ||| we describe the selection of the processes ZZ — > £ + £~£ + £~, ZZ — > £ + £~i>V, 
ZZ — > qq£ + £~ , ZZ — > qqvV, and ZZ — ► qqqq, where £ + £~ denotes a charged lepton pair of opposite 
charge and qq any of the five lightest quark-antiquark pairs. We also describe analyses of ZZ — > bb£ + £~ , 
ZZ — ► bbvV and ZZ — > qqbb which use b-tagging methods similar to those used in the OPAL Higgs 
search [||. The use of b-tagging improves the separation of the Z-pair signal from background and 
allows us to check the bb content of our Z-pair sample for consistency with the Standard Model. 
The description of the individual selections is followed by a discussion of possible systematic errors 
(Section |3.6|) . In Section |] the selected events are used to measure the Z-pair cross section. Then the 
cross section and angular distribution are compared with the Standard Model predictions and limits 
on anomalous neutral current triple gauge couplings are derived. 




e+ Z e + Z 



Figure 1: NC2 Feynman diagrams for the process e + e — * ZZ leading to a final state with four 
fermions. 
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2 Data analysis and Monte Carlo 



The OPAL detectorfl, trigger and data acquisition system are described fully elsewhere |7|-[Tl|. Our 
analyses use approximately 55 pb~ x of data collected at center-of-mass energies between 181-184 GeV 
and approximately 178 pb _1 collected at center-of-mass energies near 189 GeV. The corresponding 
luminosity- weighted mean center-of-mass energies are 182.62 ±0.05 GeV and 188.63 ±0.04 GeV [|l~2^. 
The luminosity was measured using small-angle Bhabha scattering events recorded in the silicon- 
tungsten luminometer [11, 1^, 14 1 and the theoretical calculation of Reference [ji5f| . The overall error 



on the luminosity measurement amounts to less than 0.5% and contributes negligibly to our cross- 
section measurement error. 

Selection efficiencies and backgrounds were calculated using Monte Carlo simulations. All events 



were passed through a simulation |16[ of the OPAL detector and processed as for data. We define 
the ZZ cross section as the contribution to the total four-fermion cross section from the NC2 Z-pair 
diagrams shown in Figure |j. All signal efficiencies given in this paper are with respect to these Z-pair 
processes. Contributions from all other four-fermion final states, including interference with NC2 
diagrams, are considered as background. For studies of the signal efficiency we have used grc4f [|l7j], 



YFSZZ @ and PYTHIA ||. 



Backgrounds are simulated using several different generators. PYTHIA is used to simulate two- 
fermion final states such as e + e~ — > Z*(nj) — > qq(«7) and e + e _ — > 7* (nj) — > qq(n7), where (nj) 
indicates the generation of one or more initial state photons. HERWIG (2(| and KK2f [^TJ are used 
as checks for these final states. These two-fermion generators include gluon radiation from the quarks 
which produce qqg, qqqq and qqgg final states. The grc4f generator, with the contribution exclusively 
due to NC2 diagrams removed, is used to simulate other four-fermion background. KORALW [E^] 
and EXCALIBUR are used as checks of the four-fermion background. Multiperipheral ("two- 
photon") processes with hadronic final states are simulated by combining events from PYTHIA, for 
events without electrons scattered into the detector, and HERWIG pCfl , for events with electrons 
scattered into the detector. For the qqe + e~ final state, TWOGEN is used to simulate two-photon 
events with both the electron and positron scattered into the detector. The Vermasercn [^] generator 
is used to simulate multiperipheral production of the final states e + e~£ + £~ . 

To avoid background from four-fermion final states mediated by Z7*, our selections were op- 
timized to select simulated events with masses, mi and m-2, that satisfy mi + mi > 170 GeV and 
|mi - m 2 | < 20 GeV. At 189 GeV (183 GeV) more than 90% (80%) of the events produced via the 
NC2 diagrams are contained in this mass region. Events from the NC2 diagrams dominate in this mass 
region except for final states containing electron pairs. Backgrounds in these samples from two-photon 
and electroweak Compton scattering (e7 — > eZ) processes |2(J are reduced by using electrons detected 
in the electromagnetic calorimeters with | cos# e | < 0.985, where 6 e is the polar angle of the electron. 



3 Event selection 

In the following subsections we describe event selections which exploit every decay mode of the Z 
boson. Our selections cover all ZZ final states except vvvv and t + t~vV. In hadronic final states, the 
energy and direction of the jets are determined using reconstructed tracks and calorimeter clusters 
using the correction for double counting described in Reference ||. In the qql + £~, qqqq and qqbb 
analyses four-constraint (4C) and five-constraint (5C) kinematic fits are used. The 4C fit imposes 
energy and momentum conservation. In the 5C fit the added constraint requires the masses of the two 

x OPAL uses a right-handed coordinate system in which the z axis is along the electron beam direction and the x axis 
is horizontal. The polar angle, 0, is measured with respect to the z axis and the azimuthal angle, <f), with respect to the 
x axis. 
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candidate Z bosons to be equal to one another. For final states with one or more Z bosons decaying to 
tau pairs, the energy and total momentum of the tau leptons are obtained by leaving the reconstructed 
direction of the four fermions fixed and scaling the energy and momentum of each of the fermions to 
obtain energy and momentum conservation. The scaled values of the tau momentum and energy are 
then used in the subsequent steps of the analysis. In the qqr + T~ and bbr + r _ final states, subsequent 
kinematic fits are effectively 2C and 3C fits. 

3.1 Selection of ZZ -> £+£'£+£' events 

Z-pair events decaying to final states with four charged leptons (£ + l~l^l~) produce low multiplicity 
events with a clear topological signature that is exploited to maximize the selection efficiency The 
l + £~l + £~ analysis begins by selecting low multiplicity events (less than 13 tracks or clusters) with 
visible energy of at least Q.2^fs and at least one track with momentum of 5 Ge V or more. Using a cone 
algorithm, the events are required to have exactly four cones of 15° half angle containing between 1 
and 3 tracks. Cones of opposite charge are paired^ to form Z boson candidates. 

Lepton identification is only used to classify events as background or to reduce the number of cone 
combinations considered by preventing the matching of identified electrons with identified muons. 
Electrons are identified on the basis of energy deposition in the electromagnetic calorimeter, track 
curvature and specific ionization in the tracking chambers. Muons are identified using the association 
between tracks and hits in the hadron calorimeter and muon chambers. 

To reduce background from two-photon events with a single scattered electron detected, we elim- 
inate events with forward going electrons (backward going positrons) with the cut cos# e - < 0.85 
(cos 6^+ > —0.85). Here 9 e - ( e +) is the angle of the electron (positron) with respect to the incoming 
electron beam. Background from partially reconstructed qq(n7) events and two-photon events is re- 
duced by requiring that most of the energy is not concentrated in a single cone, E v i s — E™^, > 0.2y/s. 
Here E v [ s is the total visible energy of the event and i^ne * s the energy contained in the most energetic 
cone. 

The invariant masses of the lepton pairs are calculated in three different ways which are motivated 
by the possibility of having zero, one or two tau pairs in the event. The events are classified according 
to the number of tau pairs in the event, (i) Events with E v i s > 0.9-y/s are treated as e + e~e + e~, 
e + e~ / u + / u _ or fi + /i~ p + \x~ events. We also treat all events with |cos6* m i ss | > 0.98 (6> m i ss is the polar 
angle associated with the missing momentum in the event) as e + e~e + e~, e + e~p + p~ or p + pT p + \x~ 
events to maintain efficiency for Z-pairs with initial state radiation. In these events there are no 
missing neutrinos and the mass of each pair of lepton cones is evaluated, (ii) Events failing (i) with a 
cone-pair combination that has energy exceeding 0.9mz are tried as an e + e~T + T~ or /x + p~T + r~ final 
state. The mass of the tau-pair system is calculated from the recoil mass of the presumed electron or 
muon pair, (iii) Any remaining combinations are treated as t + t~t + t~ final states. The momenta of 
the tau leptons are determined with the scaling procedure described in the introduction to Section |3] 
and the invariant masses of the cone pairs are evaluated using the scaled momenta. In any event 
with more than one valid combination, each combination is tested using the invariant mass cuts listed 
below. 

To reduce the combinatorial background, combinations with pair masses closest to mz are selected. 
In events with one or more combination satisfying \mz — mu\ > O.lmz the cone-pair combination 
with the smallest value of (mz — rnie) 2 + (mz ~ m e'£') 2 1S selected for further analysis. In the other 
combinations, the combination with the smallest value of \rnz~ mu\ or \m%— mi>(i\ is selected. The final 
event sample is then chosen with the requirement mm + m^gi > 160 GeV and \m^ — min>\ < 40 GeV. 

2 Two-track cones are assigned the charge of the most energetic track if the momentum of one track exceeds that of 
the other by a factor of 4. Events with a cone which fails this requirement are rejected. 
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The signal detection efficiency, averaged over all £ + £~£ + £~ final states is given in Table p] (line a). The 
efficiency for individual final states range from 30% for t + t~t + t~ to more than 70% for fi + fi + /j>~ . 
The invariant masses of all cone pairs passing one of the selections are shown in Figure ||a. One 
candidate is found in the 183 GeV data and one candidate is found in the 189 GeV data. 



3.2 Selection of ZZ — > t+l~vv events 

The selection of the e + e - ^77 and /U + \x~vv final states is based on the OPAL selection of W pairs 



decaying to leptons |27|. The mass and momentum of the Z boson decaying to vV are calculated using 
the beam energy constraint and the visible decay of the other Z boson to a charged lepton pair. A 
likelihood selection based on the visible and recoil masses as well as the polar angle of the leptons, is 
then used to separate signal from background. 

The e + e~vV selection starts with OPAL W-pair candidates where both charged leptons are classi- 
fied as electrons. Each event is then divided into two hemispheres using the thrust axis. The highest 
momentum charged (leading) track is selected from each hemisphere. The sum of the charges of these 
two tracks is required to be zero. The determination of the visible mass, m v i s , and the recoil mass, 
^recoib is based on the energy as measured in the electromagnetic calorimeter and the direction of the 
leading tracks. 

Three variables were chosen for the likelihood selection: QcosO, where 9 is the angle of the 
highest momentum charged track and Q is its charge, the normalized sum of visible and recoil masses 
(m v i s + m reco ii)/i/s and the difference of visible and recoil masses, m v i s — m reco n. The performance 
of the likelihood is improved with the following preselection: —25 GeV < m v ; s — m reco ii < 15 GeV and 
(m v \ s + m iecoi \) / ^ > 0.90. One event with C c + e - vV > 0.60 is selected (see Table [l] (line b) and 
Figure ||a) . 

The (j7 [x~vTi selection starts with the OPAL W-pair candidates where both charged leptons are 
classified as muons. The selection procedure is the same as for the e + e~i'T ! final states except that 
m v i s , m reco ii, and E v [ s are calculated from the the momentum of the reconstructed tracks of the Z 
boson decaying to muon pairs. The likelihood preselections —25 GeV < m v i s — m reco ii < 25 GeV and 
(m v i S + m recoi i)/ i/s > 0.90 are applied. Two events with C^+^-vv > 0-60 are selected (see Table [j] 
(line c) and Figure ||b). 



3.3 Selection of ZZ — > qq£+€ events 

The lepton pairs in the qqe + e~ and qq~/U + fj,~ final states have a distinctive signature making possible 
selections with high efficiencies and a low background contamination. In the qqr + r~ final state, the 
decay of the tau leptons produces events which are more difficult to identify. The identification of this 
final state exploits the missing momentum and missing energy carried away by the neutrinos produced 
in the decay of the tau lepton. 



3.3.1 Selection of ZZ — > qqe+e and ZZ — ► qq/x^/x events 

The selection of qqe + e~ and qqfi + n~ final states requires the visible energy of the events to be greater 
than 90 GeV and at least six reconstructed tracks. Among all tracks with momenta greater than 2 GeV, 
the highest momentum track is taken as the first lepton candidate and the second-highest momentum 
track with a charge opposite to the first candidate is taken as the second lepton candidate. Using the 
Durham 28] jet algorithm, the event, including the lepton candidates, is forced into four jets and the 
jet resolution variable that separates the three-jet topology from the four-jet topology, 2/34, is required 
to be greater than 10 -3 . Excluding the electron or muon candidates and their associated calorimeter 
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Figure 2: (a) Invariant masses of £ + £~~£ + £~ cone pairs for the data at y/s = 183 GeV and y/s = 
189 GeV. Invariant mass pairs for the (b) qqe + e~ (c) qqu + //~ and (d) qqr + T _ data. The dashed lines 
show the final invariant mass cuts. 
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Figure 3: Likelihood discriminant used at y/s = 189 GeV for (a) the e + e uv selection, (b) the 
fi + fi~vV selection, (c) the qquV selection, (d) the qqqq selection and (e) the qqbb selection. 
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clusters, the rest of the event is forced into two jets. The 4C and 5C fits to the two lepton candidates 
and the two jets are required to converge.^ 

In the qqe + e~ selection no explicit electron identification is used. Electron candidates are selected 
by requiring the sum of the electromagnetic cluster energies E\ + Ei associated to the electrons to be 
greater than 70 GeV and the momentum of the most energetic electron track to exceed 20 GeV. We 
also reject the event if the angle between either electron candidate and any other track is less than 5°. 

In the qq/i + ^~ selection the muons are identified using (i) tracks which match a reconstructed 
segment in the muon chambers, (ii) tracks which are associated to hits in the hadron calorimeter or 



muon chambers [27], or (iii) isolated tracks associated to electromagnetic clusters with reconstructed 
energy less than 2 GeV. No isolation requirement is imposed on events with both muon tracks passing 
(i) or (ii). Events with at least one muon identified with criterion (iii) are accepted if both muon 
candidates in the event have an angle of at least 10° to the nearest track. We require the sum of the 
momenta of the two leptons to be greater than 70 GeV. 

Z-pair events are separated from Z7* background by requiring the fitted mass of the 5C fit to 
be larger than 85 GeV and the invariant masses my and m qq obtained from the 4C fit to satisfy 
(mu + m qq ) > 170 GeV and — m qq | < 30 GeV. Figure^ (^c) shows the distribution of m ee {m^) 
and m qq before the cuts on the masses from the 4C and 5C fits. 

After all cuts the selection efficiency^ for qqe + e~ signal events is (55.2 ± 2.7)% at 183 GeV and 
(65.1 ±2. 9)% at 189 GeV. The errors on these efficiencies include the systematic errors (see Section |3~6| ) . 
No candidate is observed at 183 GeV. Six candidate events are found after all cuts in the data taken 
at 189 GeV. In Table [l] (lines d and e) we give the efficiency, background and observed number of 
events. The largest source of background after all cuts is from Z7* mediated qqe + e~ events and from 
two-photon events. 

For the qqfj, + /j,~ final state the selection efficiency is (66.9 ± 2.6)% at 183 GeV and (72.3 ± 2.8)% 
at 189 GeV. Three events are observed in the 189 GeV data sample and none in the 183 GeV data 
sample. The errors on these efficiencies include the systematic errors given below. In Table |l| (lines f 
and g) we give the efficiency, backgrounds and observed number of events. The background after all 
cuts is expected to come from e + e~ — > Z7* — > qq// + ^~ events. 



3.3.2 Selection of e+e — ► ZZ — ► qqr+r events 



The qqr + T - final state is selected from a sample of events with track multiplicity greater or equal 
to six. Events which have been selected as e + e~ — > ZZ — ► qqe + e~ or e + e~ — > ZZ —* qq/i + /U _ are 
excluded from this selection. The tau-lepton candidates are selected using an artificial neural network 
algorithm which is described in detail in Reference |]]. The tau candidate with the highest neural 
network output value is taken as the first candidate. The second best candidate is required to have its 
charge opposite to the first candidate and the highest output value among all remaining candidates. 
If a second candidate cannot be found the event is rejected. 

Motivated by the presence of neutrinos in the final state, the visible energy of the event, E v [ s , 
is required to exceed 90 GeV and the missing energy — E v i s is required to exceed 15 GeV. In 
addition, the sum of the momenta of the leading tracks from the tau-lepton decays is required to be 
less than 70 GeV. Since the direction of the missing momentum in signal events will tend to be along 
the direction of one of the decaying tau leptons, the angle a r ,miss between the missing momentum and 
a tau-lepton candidate is required to satisfy a Tjm i ss < 90° for at least one of the two tau candidates. 

3 In the context of this paper, convergence is defined as a fit probability greater than 10 -10 . 

4 Small amounts of feedthrough from other ZZ final states, in this case qqr + r~, are counted as signal. 
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The two hadronic jets are selected in the same way as in the e + e~ — > ZZ — > qqe + e~ selection. 
The initial estimate of the energy and the momenta of the tau candidates is found from the sum of 
the tracks associated to the tau by the neural network algorithm and all unassociated electromagnetic 
clusters in a cone with a half angle of 10° around the leading track from the tau decay. A 2C kinematic 
fit that imposes energy and momentum conservation (see the introduction to Section ||) is required to 
converge. A 3C kinematic fit, with the additional constraint of the equality of the fermion pair masses 
is also required to converge. 

Using the network output || for each tau lepton, a probability is calculated taking into account 
the different branching ratios, sensitivities, efficiencies and background levels for 1-prong and 3-prong 
tau-lepton decays. In the following, we combine the probabilities V\ and V2 to form a likelihood using 



v 1 v 2 + (i-r l )(i-v 2 )' KJ 

The likelihood associated with probabilities of the two tau candidates is required to satisfy C TT > 0.977. 
In addition, the common mass of the 3C fit is required to exceed 85 GeV. Using the 2C fit masses 
of the tau pair, m TT , and the quark pair, m qq , as obtained from the kinematic fit, we also require 
m-qq + Tn TT > 170 GeV and \m m — m TT \ < 60 GeV. 

After all cuts the selection efficiency for signal events is found to be (22.2 ± 1.7)% at 183 GeV and 
(26.9 ± 2.0)% at 189 GeV. One candidate event is found in the data at 189 GeV while no candidate 
is selected at 183 GeV. Figure ||d shows the masses of the candidate events before the invariant mass 
cuts. In Table |l| (lines h and i) we give the efficiencies, backgrounds and observed number of events. 



3.3.3 Selection of e+e~ -»• ZZ -»■ hh£+£~ 



Events with bb final states are selected using the algorithm described in Reference ||. The b prob- 
abilities of the two hadronic jets are combined to form a likelihood, £bb> according to Equation [l|. 
Because the qqe + e~ and qq/x + /i~ selections are pure, a relatively loose cut of £t>b > 0.2 is used to 
select the bbe + e~ and bb/x + /i~ samples. For the selections with electron and muon pairs there are 
two classes of events since the selected bb£ + £~ events are a subset of the qq£ + £~ events. In Table [l] 
(lines e and g) we give the efficiencies of the b-tagged samples with respect to the expected fraction of 
bb events. The efficiencies for samples without b-tags are given with respect to the hadronic decays 
without bb final states. 

In the bbe + e~ selection one candidate is found in the data. The selection efficiency is found to 
be (47±3)% at 183 GeV and (46±3)% at 189 GeV. Qln the bb/i + /i~ selection we find no candidate at 
183 GeV and one candidate at 189 GeV. The selection efficiency is found to be (50±3)% at 183 GeV 
and (54±3)% at 189 GeV. 

For the bbr + r~ selection the C TT cut of the qqr + r~ selection is loosened and combined with £t>b 
as follows. C TT and £bb are both required to be greater than 0.1. The bbr + r~ probability for the 
event, £bbi-r> is calculated from Equation [j] with C TT and £bb as inputs and required to exceed 0.95. 
After the cut on £bbTi-> the remaining cuts of the qqr + r _ selection are applied, giving a selection 
efficiency of (21 ± 3)% at 183 GeV and (24 ± 3)% at 189 GeV. No candidate event is found in the 
183 GeV or 189 GeV data. 

We also use an alternative jet-based bbr + T~ selection || and accept any event which passes 
either bbT + r~ analysis, but events previously selected by another qq£ + £~ selection are rejected. The 
alternative analysis uses a different approach to reconstruct the tau leptons. This event selection 

5 The efficiencies given in this section do not include feedthrough from other qq£ + l~ final states. The efficiencies 
given in Table ^ include this feedthrough. 
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consists of a set of preselection cuts and a subsequent multivariate likelihood selection. Events are 
reconstructed as four jets using the Durham algorithm. Tau-lepton candidates are sought in the four 
jets using a likelihood technique to separate real tau leptons and fakes in quark jets. The tau and b-tag 
likelihood values are combined in a bbr + T~ likelihood which is maximized to choose the b jets and 
tau leptons of the event. This bbr + T~ likelihood uses tau and b-tag likelihoods and some topological 
variables as input. Events are accepted if their likelihood exceeds 0.6. In addition, the fitted 2C 
masses are required to satisfy m qq + m TT > 170 GeV and |m qq — m TT \ < 60 GeV. 

At 189 GeV data the alternative selection has an efficiency of (30 ± 3)% for bbr + r~ events. After 
combining the two selections, the efficiency for bbr + r~ events for all cuts except those rejecting events 
found by the other qq^ + £~ selections is 40%. A similar improvement is realized for the 183 GeV data. 
The efficiency and backgrounds after rejecting events found by the other selections are given in Table [l] 
(line j). One exclusive event is selected by the combined bbr + r~ analysis at 189 GeV. 



3.4 Selection of ZZ — > qqi/i/ events 

The qqvV selection is based on the reconstruction of the Z boson decaying to qq which produces 
somewhat back-to-back jets. The selection uses contained events with a two-jet topology. The beam 
energy constraint is then used to determine the mass of the Z boson decaying to vV. The properties of 
the qq decay and the inferred mass of the vV decay are then used in a likelihood analysis to separate 
signal from background. 

Two-jet events are selected by dividing each event into two hemispheres using the plane perpen- 
dicular to the thrust axis. The number of charged tracks in each hemisphere is required to be four 
or more. The polar angles of the energy-momentum vector associated with each hemisphere, #hemii 
and #hemi2; are used to calculate the quantity cos#h = ^(cos#hcmii — cos#hemi2)- Contained events 
are selected by requiring |cos#h| < 0.80. The total energy in the forward detectors and in the for- 
ward region of the electromagnetic calorimeter (| cos#| > 0.95) is required to be less than 3 GeV. W 
decays identified by the OPAL W-pair selection are rejected; the likelihood for e + e~ — ► qq£z^ from 



Reference [27], £ww> is required to be 0.5 or less. 

An important background to our selection is qq(«7) events with photons which escape detection. 
We discriminate against these events by looking for a significant amount of missing transverse mo- 
mentum, p t . In each event, pt can be resolved into two components, pu, perpendicular to both the 
thrust axis and the beam axis and ptj, along the thrust axis and perpendicular to the beam axis. 
Since pa is based primarily on angular measurements, it is better measured than p t j. We approxi- 
mate pa as pa = \E\> sin^sin^h- Here E\> = \fs/2 is the beam energy, <f> is the acoplanarity of the 
momentum vectors of the two hemispheres and sin 9^ = \/l — cos 2 6^. The resolution on pa , a Pti , was 
parameterized as a function of thrust and cos #h using data taken at the Z resonance. The variable 
Rpu = (Pti ~ Pti)/ a pti i s use d as an input to the likelihood. Here p^ corresponds to the transverse 
momentum carried by a photon with half the beam energy which just misses the inner edge of our 
acceptance ( = i?bsin(32mrad)/2). We also use the variable cos6> m i ss , the direction of the missing 
momentum in the event, to discriminate against the qq(n7) events. 

In the final selection of events, we use a likelihood based on the following five variables: (i) the 
normalized sum of visible and recoil masses (m v ; s + m reC oii ( n ) the difference of visible and recoil 

masses (m v ; s — m reC oil), (iii) log(?/23)> where 1/23 is the jet resolution parameter that separates the two- 
jet topology from the three-jet topology as calculated from the Durham jet algorithm, (iv) cos# m i ss 
and (v) Rp H - The mass variables are useful for reducing background from W-pair production and Weis 
final states. The jet resolution parameter is useful in reducing the remaining qq£i/ final states. To 
improve the performance of the likelihood analysis we use only events with : |m v i s — m reco ;i| < 50 GeV, 
(m v i s + mrecoiO/v^ > 0.89 and R Pti > 1.2. Events are then selected using C^-p > 0.5, where £ qqj /j7 
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is the likelihood for the qqz/77 selection. The likelihood distribution of data and Monte Carlo is shown 
in Figure ||c. For the bbvV selection we require, in addition, the b-tag variable of Reference || to be 
greater than 0.65. 

The efficiencies for the qquV selection alone are (30.5 ± 2.0)% at 183 GeV and (33.9 ± 2.3)% at 
189 GeV. The errors on these efficiencies include the systematic errors discussed below in Section 
The efficiencies after considering the results of the b-tagging, as well as the number of events selected 
at the two energies are given in Table [l] (lines k and 1) . 



3.5 Selection of ZZ — ► qqqq events 



The fully hadronic channel of the Z-pair decay has the largest branching fraction of all channels (about 
50%), but suffers from large background from hadronic W-pair decays. We apply two different event 
selections, one of which is based mainly on reconstructed mass information in order to accept all 
hadronic Z-pair decays without flavor requirement, while a second analysis applies a flavor tag in 
order to select final states involving b quarks, allowing for looser requirements on the reconstructed 
boson mass. 

For both subsamples, hadronic Z*/7* — ► qq events are an important background. We therefore 
start with a common preselection based on event shape variables which is mainly aimed at reducing 



this background. We use a likelihood method, described below in Section 3.5.1, in order to choose the 
most likely jet pairing for Z-pair decays. 



3.5.1 Preselection and jet pairing 

The event selection starts from the inclusive multihadron selection described in Reference [O]. The 
radiative process e + e~ — ► Z7 — > qqy is suppressed by requiring the effective center-of-mass energy 
after initial state radiation, v?, to be larger than 150 GeV. yfs' is obtained from a kinematic fit [O] 
that allows for one or more radiative photons in the detector or along the beam pipe. The final state 
particles are then grouped into jets using the Durham algorithm |ffif| . A four-jet sample is formed by 
requiring the jet resolution parameter 2/34 to be at least 0.003 and each jet to contain at least two 
charged tracks. In order to suppress Z*/7* — > qq background, the event shape parameter C par [pOfl , 
which is large for spherical events, is required to be greater than 0.25. A 4C kinematic fit using energy 
and momentum conservation is required to converge. A 5C kinematic fit which forces the two jet pairs 
to have the same mass is applied in turn to all three possible combinations of the four jets. This fit is 
required to converge for at least one combination. The efficiencies of these preselection cuts are (86.4 
± 0.5) % and (88.9 ± 0.5) % for signal events at 183 GeV and 189 GeV, respectively. 

In order to determine which pair of jets comes from each Z, we calculate a likelihood function using 
the mass obtained from the 5C fit, the corresponding fit probability, and the difference between the 
two di-jet masses obtained from the 4C fit. In the YFSZZ simulation of Z-pair decays the correct jet 
pairing has the highest likelihood output in (86.8 ± 0.5) % of the events. This fraction rises to (93.8 
± 0.5) % for the events after the final selection. 



3.5.2 Likelihood for the inclusive ZZ — > qqqq event selection 

We use a likelihood selection with eight input variables for the selection of ZZ —* qqqq events. The 
first variable is the jet pairing likelihood described above. Excluding the jet pairing with the largest 
difference between the two di-jet masses as obtained from the 4C fit, we identify among the remaining 
two possible pairings the one for which the 5C-fit mass is closer to the W mass. We use the difference 
between this 5C-fit mass and the W mass in order to discriminate against hadronic W-pair events. 
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Two variables that are sensitive to unobserved particles along the beam direction are the fitted center- 
of-mass energy and the sum of the cosines of the polar angles of the four jets. In order to discriminate 
against Z*/7* events, we use the difference between the largest and smallest jet energies after the 
4C fit, and the angular variable j ang = E^(l — cos 9\2 cos #13 cos 023)/\/s, where E± is the smallest of 
the four jet energies, and the 6ij are the opening angles between jets i and j, with the jets ordered 
by energy. Finally we calculate from the momenta configuration of the four jets the effective matrix 
element for the QCD processes Z*/7* — > qqgg and Z*/7* — > qqqq as defined in Reference p9|| , and 



the matrix element for the process WW — > qqqq from Reference [23]. 

The distribution of the likelihood function calculated from these eight variables is shown in Fig- 
ure ||d. In order to maximize the significance of the measured cross section, assuming a 10% relative 
systematic error on the background, we place a cut on the likelihood at 0.65. This cut leads to an 
efficiency of (33.0 db 1.8)% (189 GeV) relative to all fully hadronic NC2 Z-pair final states. In the data, 
52 events are selected. We expect a total of 57.0 events from Standard Model processes, of which 17.6 
originate from the Z-pair signal, while 27.0 events are expected from hadronic W-pair decays and 12.4 
events from hadronic two-fermion processes. At 183 GeV the selection efficiency for fully hadronic 
Z-pair decays is (25.3 ± 1.3)%. In the data 8 events are selected. The Standard Model expectation is 
1.7 signal events and 7.2 background events. 



3.5.3 Likelihood for ZZ — ► qqbb event selection 

Jets originating from b-quarks are selected using the same b-tagging algorithm used in the bb£ + £~ 
and bbz^I7 selections. We evaluate the probability for each of the four jets to originate from a primary 
b quark, and use the two highest probabilities as input variables for a likelihood to select ZZ — ► qqbb 
events. In addition, we use the parameters 2/34, C par , the difference between the largest and smallest 
jet energies and the output of the jet pairing likelihood. We also use the fit probabilities of a 5C 
kinematic fit which constrains one boson mass to the Z mass, and the probability of a 6C fit which 
forces both masses to be equal to the W mass. 

Figure ||e shows the distribution of the likelihood function calculated from these eight variables for 
the preselected events. The signal likelihood is required to be larger than 0.80 for both 183 GeV and 
189 GeV data. After the likelihood selection, we perform an additional cut on the mass obtained from 
the 5C-mass fit in the most likely jet pairing, which is required to be larger than 86 GeV. This choice 
of the cuts on the likelihood and on the 5C-fit mass was made by maximizing the expected statistical 
significance of signal over background. The final efficiency is (35.1 ± 2.4)% for the 183 GeV data and 
(38.6 ± 2.7)% for the 189 GeV data. The observed number of events, expected signal and background 
are given in Table [l] (line o) . 



3.5.4 Combination 

To account for overlap between the qqqq and qqbb selections we divide the data into three logical 
classes, exclusive qqqq events (Table [j] (line m)), exclusive qqbb (line n) events and the overlapping 
region (Table |l| (line o)). In Table |l| (line m) we give the efficiency relative to all fully hadronic final 
states without b-quarks. The other efficiencies are given relative to fully hadronic final states with 
b-quarks. 



3.6 Selection systematic errors 

Systematic errors have only a modest effect on our final result because of the large statistical error 
associated with the small Z-pair cross section. 
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Detector effects can best be studied by comparing calibration data taken at the Z resonance with 
a simulation of the same process. These comparisons are important for final states such as £ + £~vV, 
qqvV, and qqqq, where the tight cuts are needed to separate signal and background. In these cases, 
we add additional smearing to the energy and momentum of the simulated events to match data and 
simulation. We then apply the same smearing to the signal and background Monte Carlos and then 
correct our efficiencies and background accordingly. The full difference is used as the systematic error 
in these cases. At y/s = 189 GeV, these differences give relative systematic errors on the efficiency of 
2.5% for the e + e~^I7 final state, 5.2% for the /i + /x~z^I7 final state and 3.8% for qquu final state. In the 
qqqq inclusive analysis, where uncertainties on the reconstructed angles are also important, similar 
studies lead to a relative detector systematic error of 6%. 

Detector systematic errors for the qq£ + £~ and £ + £~£ + £~ selections without r-pairs in the final 
state are small because of the good separation of signal and background. In the £ + £~£ + £~~ final state 
the largest effect (3%) is from modeling of the multiplicity requirement which is important for final 
states containing r-pairs. In the qqr + r _ selection, the systematic uncertainties in the efficiencies were 
determined by overlaying hadronic and tau decays taken from Z resonance data giving a contribution 
to the systematic error of 6.0%. When the final states are combined to determine the cross section 
and limits on the anomalous triple gauge couplings, we assume a common relative systematic error of 
3% on the efficiencies. 

Another important detector effect comes from the simulation of the variables used by the OPAL 
b-tag which is discussed in Reference [||. We allow for a common 5% error in the efficiency of the 
b-tag, consistent with our studies on Z resonance data and Monte Carlo. 

These detector effects were propagated through to our background errors. 

In each channel the signal and background Monte Carlo generators have been compared against 
alternative generators. In almost all cases the observed differences are consistent within the finite 
Monte Carlo statistics and the systematic error has been assigned accordingly. One notable exception 
is in the qqbb background where differences in the PYTHIA, grc4f and EXCALIBUR simulation of 
the W-pair background are as large as 20%. In this case the full difference has been assigned as the 
systematic error. 

In the £ + £~£ + £~ channels, which has a large background from two-photon events, we have com- 
pared the number of selected events at an early stage of the analysis with the Monte Carlo prediction 
and based our background systematic error on the level of agreement. This results in 20% systematic 
error. 

4 Results 

We combine the information from all of the analyses reported above using a maximum likelihood fit 
to determine the production cross section for e + e~ — ► ZZ. The information which was used in the fit, 
as well as the Standard Model prediction for Z-pair production, is summarized in Table |l|. For each 
channel the table gives the number of events observed, n b S) the Standard Model prediction for all 
events, nsM, the expected signal, nzz, the expected background, nb ac k> the efficiency e c h a n> an d the 
integrated luminosity, L m t. Bzz is the branching ratio of Z-pairs to the given final state, calculated 
from Z resonance data 0. In the table we give the overlap between the b-tag and non-b-tag analyses. 
Possible overlap between qqqq and qq£ + £~ has been studied, and found to be an order of magnitude 
smaller than the overlap of qqqq and qqbb and has therefore been ignored. 

The cross section at each energy is determined with a maximum likelihood fit using a Poisson 
probability density convolved with Gaussians to describe the uncertainties on efficiencies and back- 
grounds. The expected number of events in each channel, fj, e , as function of the Z-pair cross section, 
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0.30 ± 0.02 


0.219 


54.5 


1 


qqz/77 & bbz^I7 


1 


0.42 ± 0.06 


0.25 ± 0.02 


0.17 ± 0.06 


0.31 ± 0.03 


0.061 


54.5 


m 


nTTnTT a/ n?7l~i}"> 
mu vj u 


8 


7 43 ± 70 


1 03 ± 07 


6 39 ± 69 

\J . <J i > _1_ VJ . VJ sJ 


26 ± 02 

KJ • 4-/KJ 1 KJ • KJ £-1 


0.300 


54.7 


n 
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1 co i n oi 

l.OU _l_ VJ.tJJ- 
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I.I I vJ.tJ-L 


o 1 b + n oi 

\J.±KJ _l_ vJ.VJ-L 
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VJ . ' J ' J 


54 7 


KJ 
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±. L ±± _1_ VJ.XiJ 
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VJ.OVJ _l_ VJ.-LiJ 
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kj.£kj _1_ VJ.VJi 


935 

VJ . LsKJKJ 
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\J • KJ KJ — 1 — KJ • KJ £J 


35 + 21 

V 7 . ' J " 7 — 1 — V.' . — .1. 


58 zb 02 

KJ iKJKJ — 1 — U • U 


0.010 


182.1 


b 




1 


1.30 ± 0.17 


0.72 ± 0.09 


0.58 ± 0.15 


0.46 ± 0.06 


0.013 


181.5 




1I~ uv 

LA/ LA/ IS IS 


2 


97 ± 15 


54 ± 08 


0.44 ± 0.13 


35 ± 05 

V 7 . ' ' KJ — 1 — V/ • V / KJ 


0.013 


181.5 


d 




5 


3 51 + 38 


2 79 ± 13 

— ■ . 1 — 1 — VJ • X KJ 


72 ± 36 

KJ • 1 £1 — 1 — VJ • KJ \J 


69 ± 03 

KJ • VJ tJ — 1 — KJ • KJ KJ 


0.037 


174.7 


e 


qqe^e - & bbe^e - 


1 


0.68 ± 0.07 


0.59 ± 0.05 


0.09 ± 0.06 


0.52 ± 0.05 


0.010 


174.7 


f 


VJ VJ i_J/ X-*/ tJO PL' J-J X-*/ J--*/ 


2 


3 22 ± 12 


3 05 ± 10 

' J . V 7 " ./ — 1 — VJ • ±. KJ 
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KJ • VJ — 1 — VJ • VJ VJ 


75 ± 02 

KJ • 1 KJ — 1 — KJ • KJ £J 


0.037 
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& 




1 


0.73 ± 0.04 


0.71 ± 0.04 


0.02 ± 0.02 


0.63 ± 0.05 


0.010 


174.7 
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i 
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1 OO _|_ Q 1 O 
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1 74 7 
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qqT^T - Sz bbT~^T~ 





30 ± 06 


25 ± 04 

\J * w — I — V / • V / -i- 


04 ± 04 


0.22 ± 0.04 


0.010 


174.7 


j 


bbr + r _ & qqr+T~ 


1 


0.29 ±0.05 


0.27 ±0.04 


0.02 ±0.01 


0.24 ±0.04 


0.010 


174.7 


k 


qqi/I7 & hhuu 


20 


16.1 ± 1.2 


8.5 ±0.7 


7.51 ±0.94 


0.35 ±0.03 


0.219 


174.2 


i 


qquV &: hhvu 


5 


2.18 ±0.24 


1.90 ±0.20 


0.28 ±0.14 


0.29 ±0.03 


0.061 


174.2 


m 


qqqq & qqbb 


45 


50.1 ±4.0 


12.5 ±0.8 


37.6 ±3.9 


0.38 ±0.03 


0.300 


174.7 


n 


qqbb & qqqq 


11 


6.74 ± 0.75 


3.57 ±0.23 


3.16 ±0.72 


0.14 ±0.01 


0.235 


174.7 


o 


qqbb & qqqq 


7 


7.87 ±0.64 


5.49 ±0.32 


2.38 ±0.56 


0.21 ±0.02 


0.235 


174.7 



Table 1: Observed number of events, n b s , the total Standard Model expectation, nsM> the expected 
number of Z-pairs, nzz, background expectation, nb ac k> and efficiencies, e c h an ) for the 183 GeV and 
189 GeV data. Bzz is the product branching ratio for the final state which is calculated directly from 
Z resonance data ||. Note that the efficiencies for selections with b-tags are given relative to the 
fraction of hadronic final states which contain a Z boson decaying to bb. For selections of events with 
hadronic final states, but without b-tags, the efficiencies are relative to those hadronic final states 
which do not include a Z boson decaying to bb. The errors include contributions from the common 
systematic errors. 
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" u z 180 185 190 195 200 205 

Energy (GeV) 

Figure 4: The OPAL measurements of the NC2 Z-pair production cross section. The shaded band 
shows the prediction, using the default settings, of the YFSZZ Monte Carlo for the total cross section. 
The band indicates a theoretical scale uncertainty of ±3%. 
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ozz, is given by 

(2) 

The efficiencies, e c han, include the effects of off-shell Z bosons that are produced outside of our kine- 
matic acceptance. Our main result, the NC2 Z-pair cross sections obtained from the fits, is 

<rzz(183 GeV) = 0.12 Iffistat.) i^l(syst.) pb 

azz(189 GeV) = 0.80 t{ji|(stat.) ±g;g|(syHt.) pb. 

The 183 GeV result corresponds to a 95% C.L. upper limit (normalized to the region <tzz > 0) of 
0.55 pb. The comparison of these measurements with the YFSZZ prediction, using the coupling of 
Z bosons to electrons measured at the Z resonance, is shown in Figure (§. The results are consistent 



with the YFSZZ prediction and with the measurements presented in Reference [31]. 



Our results assume that, apart from the backgrounds discussed above, only ZZ production con- 
tributes inside our kinematic region. Possible effects of Higgs boson production have been ignored. 

In order to check the results for consistency with the expected fraction of bb final states, we 
perform a second fit where the branching ratio of the Z boson to bb is a free parameter. The relative 
branching ratios of the Z to other fermion pairs are fixed to their measured values. The resulting 
branching ratio and cross section for the 189 GeV data are Br(Z — > bb) = 0.21^QQg(stat.) ± 0.01(syst.) 
and <7zz = 0.75^0^4 (stat.) ^{^(syst^pb. The measured branching ratio can be compared with the 
world average as measured at the Z resonance of Br(Z ^ bb) = 0.1516 ±0.0009 §§. The cross section 
measured without constraining the branching ratio is also consistent with our main result and the 
YFSZZ prediction. The smaller error on the main result cross section illustrates the advantage of 
classifying the hadronic systems as bb or non-bb. The 183 GeV data sample is too small to extract a 
meaningful value of the Z — ► bb branching ratio. 

Limits on anomalous triple gauge couplings were set using the total cross section and the | cos#z| 
distribution of our data. Here 6z is the polar angle of the Z bosons produced. In this study we varied 
the real and imaginary parts of the ZZZ and ZZ7 anomalous couplings parameterized by the form 
factors /4 ZZ , / zzz , and /^ Z7 as defined in Reference || and implemented in the YFSZZ Monte 
Carlo. The real and imaginary parts of each coupling were varied separately with all others fixed to 
zero. 

For this study we consider the effect of the anomalous couplings on the total cross section at 
yfs = 183 GeV and on the cross section in four bins of |cos#z| at = 189 GeV. The selection 
efficiencies for all final states are parameterized as function of anomalous couplings. At ^/s = 189 GeV 
the parameterization is done separately for each bin in | cos #z I • For values of the anomalous couplings 
larger than unity, much of the production of the final state fermions occurs at | cos#| ~ 1 where the 
efficiency for most channels is reduced by a factor of ~ 0.5. An uncertainty of 10%, dominated by 
Monte Carlo statistical errors, is assigned to the correction we apply to these efficiencies. 

The 95% C.L. limits on the anomalous couplings obtained from the maximum likelihood fit are 
given in Table |2[ With the exception of the i?e{/ zzz } coupling, the limits are insensitive to the sign 
and complex phase of the couplings. 



5 Conclusion 

The production cross section of e + e~ — ► ZZ has been measured using the final states £ + £~£ + £~, 
£ + £~vl7, qc[£ + £~ , qqvu, and qqqq. The number of observed events, the background expectation from 
Monte Carlo and the calculated efficiencies have been combined to measure the production cross 
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Coupling 


95% C.L. Lower Limit 


95% C.L. Upper Limit 


Re{f! ZZ } 


-2.1 


2.1 


Im{ff zz } 


-2.1 


2.1 


Re{fP Z } 


-6.2 


4.4 


Im{fF Z } 


-6.4 


6.4 


M/f 7 } 


-1.2 


1.2 




-1.2 


1.2 


Re{f ZZ "} 


-3.9 


3.6 


Im{f ZZl } 


-3.8 


3.9 



Table 2: The 95% confidence level limits on possible anomalous triple gauge couplings. 

section of the process e + e~ — > ZZ. Our measured cross sections include the effects of background and 
efficiency uncertainties. 

We have determined the cross section for e + e~ — > ZZ separately at average center-of-mass energies 
of 182.62 ± 0.05 GeV and 188.63 ± 0.04 GeV. The NC2 Z-pair cross sections were determined to be 

ctzz(183 GeV) = 0.12 ±g;^(stat.) ^(syst.) pb 

azz(189 GeV) = 0.80 Igif (frtat.) ^(syst.) pb. 

At the lower center-of-mass energy, the 95% C.L. upper limit on the cross section is 0.55 pb. The 
measurements at both energies are consistent with the Standard Model expectations. No evidence is 
found for anomalous neutral current triple gauge couplings. The 95% confidence level limits are listed 
m Table |. 
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